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要 :为 了 研究 矩形 腔 体 侧 向 正弦 周期 加 热 条 件 下 对 流 扰 动 的 成 长 和 传 热 特 性 ,本 研究 对 流体 力 


学 方程 组 进行 了 数值 模拟 。 结 果 表 明 : 随 着 格拉 晓 夫 数 G7 的 增加 ,对 流 扰 动 的 最 大 振幅 4,,, 的 成 
长 率 变 大 ,线性 成 长 阶段 的 时 间 变 短 。 对 于 高 宽 比 4 = 10 , 腔 体 宽度 d = 2 cm 的 腔 体 , 普 朗 特 数 


Pr 


=6. 949 的 流体 ,对 流 扰动 的 成 长 率 y, 随 着 格拉 晓 夫 数 Gr 交 化 的 经 验 式 为 y, =9 x10 Gr ; 


对 于 高 宽 比 4 = 10 、 腔 体 宽度 d = 6 cm 的 腔 体 , 普 朗 特 数 Pr = 0.703 的 流体 ,成 长 率 yn 随 着 格拉 
晓 夫 数 Gr 的 变化 关系 式 为 y= 8 x10G6r"? 。 当 Pr = 0.0272 时 , 热 壁面 平均 努 塞 尔 数 Nu 和 格 
拉 晓 夫 数 Gr 的 关联 式 为 Nu = 0.00002Gr*"”; 当 Pr = 6.949 时 , Nu = 0.00024GCrr 。 随 着 格拉 了 晓 
夫 数 Cr 数 增加 , 右 壁 面 的 传 热能 力 增强 。 

关键 词 :矩形 腔 体 ; 侧 向 正弦 周期 加 热 ; 对 流 扰 动 ;成 长 率 ; 传 热能 力 
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Abstract :In order to study the growth of convective disturbance and heat transfer characteristics of 
rectangular cavity with sinusoidal periodic heating, the hydrodynamic equations were numerically 


simulated. The results showed that with the increase of Grashof number Gr , the growth rate of maximum 


amplitude A,,. of convective disturbance increases while the time of linear growth phase decreases. For a 


cavity with aspect ratio A = 10 and a widthd = 2 cm and for the fluid with Prandtl number Pr = 6. 949 ， 
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the empirical formula of the growth rate y,, of convective disturbance varying with Craschof number is y。= 
9 x10 6r .Foracavity with4 = 10 andd =6cmandforthe fluid with Pr = 0.703 , the empirical 
formula of the growth rate yuvaring with Graschof number isy，= 8 x10 6. AtPr = 0.0272 , the 


correlation between the _ average Nusselt number Nu of the hot wall and Graschof number is Nu = 
0. 000 02Gr™™ ; at Pr = 6.949 ,the correlation is Nu = 0.00024Gr'™ .The heat transfer capacity of the 


right wall increases with the increase of Graschof number. 


Key words:rectangular cavity ;sinusoidal periodic heating from side; convective disturbance; growth rate; 


heat transfer capacity 


从 底部 加 热 的 矩形 腔 体内 的 Rayleigh-Bénard 对 
流 已 经 有 一 百 多 年 的 历史 '"“ 。 对 它 的 研究 由 最 初 
的 实验 理论 分 析 发 展 到 现在 的 基于 流体 力学 方程 
组 的 数值 模拟 “” ,一 直 吸引 着 大 量 的 研究 者 ,并 已 
经 获得 了 很 大 的 进展 。 特 别 在 混合 流体 分 离 比 为 负 
值 情况 下 发 现 了 行 波 现象 ,导致 了 对 该 问题 的 广泛 
深入 研究 。 不 同 结构 的 局 部 行 波 *”、 扎 动 行 
波 "具有 缺陷 的 行 波 汪 、 双 局 部 行 波 的 对 流 斑 
图 "局 部 行 波 的 局 部 结构 与 周期 性 等 一 系列 
现象 已 被 发 现 。 在 均匀 加 热 的 基础 上 ,为 适用 实际 
应 用 中 热源 不 均匀 波动 的 特点 ,科研 人 员 进 一 步 考 
虑 了 底部 局 部 加 热 和 底部 周期 加 热 对 对 流 斑 图 的 影 
响 , 以 及 水 平流 动 对 对 流 结 构 及 斑 图 形成 的 影 
响 "3 沾 。 这 些 研究 都 是 基于 Rayleigh-Bénard 对 流 的 
基本 模型 ,揭示 了 非 平 衡 对 流 系 统 的 复杂 而 丰富 的 
动力 学 特性 。 

对 于 和 矩 形 腔 体 ,如 果 从 一 侧 加 热 ,这 将 形成 另外 
一 种 自然 对 流 系统 。 文 献 [14 ] 最 早 考 虑 了 这 样 一 
个 系统 。 文 献 [15 ] 对 侧 向 加 热 的 对 流 系统 进行 了 
实验 研究 。 文 献 116-20 ] 研 究 了 侧 向 加 热 的 对 流 系 
统 中 的 非 恒定 自然 对 流传 热 特 性 温度 振动 .垂直 
边界 层 中 的 对 流 不 稳定 性 和 双 层 结构 。 文 献 [21] 
进一步 探讨 了 侧 向 加 热 的 对 流 系统 中 的 动力 学 特 
性 。 文 献 [22 ] 回顾 了 侧 加 热 腔 内 的 自然 对 流 的 研 
究 , 综 述 了 其 进展 。 对 侧 向 加 热 的 对 流 系统 的 研究 
既 有 理论 价值 ,同时 又 具有 实际 应 用 意义 。 大 部 分 
对 侧 向 加 热 的 对 流 系统 的 研究 都 是 基于 小 高 宽 比 
侧 均 名 加 热 腔 体 。 为 了 探讨 侧 向 局 部 周期 加 热 及 侧 
壁面 正弦 周期 加 热 对 腔 体 内 对 流 结 构 的 影响 ,文献 
[23-28] 进 行 了 有 益 的 探索 。 发现 加 热 方 式 明显 影 
响 对 流 斑 图 的 结构 和 对 流 斑 图 的 形成 。 文 献 129 | 
对 倾斜 腔 体 对 流 进行 了 研究 ,并 发 现 了 一 种 对 流 斑 
图 , 它 有 着 许多 不 同 的 动力 学 特性 。 有 必要 研究 较 


大 高 宽 比 腔 体 侧 壁 面 正弦 周期 加 热 对 对 流 扰动 的 成 
长 和 传 热 特性 的 影响 。 

本 研究 通过 对 流体 力学 方程 组 的 数值 模拟 , 研 
究 队 形 舱 体 侧 向 正弦 周期 加 热 条 件 下 对 流 扰 动 的 
成 长 和 传 热 特性 。 针 对 高 宽 比 4 = 10 , 腔 体 宽度 
d =2 cm 的 腔 体 , 普 朗 特 数 Pr = 6. 949 的 流体 ,和 高 
宽 比 4 = 10 、 腔 体 宽度 d = 6 cm 的 腔 体 , 普 朗 特 数 
Pr = 0.703 的 流体 ,探讨 其 对 流 扰动 的 成 长 率 y; 随 
着 格拉 晓 夫 数 Gr 变化 的 经 验 式 。 发 现 随 着 格拉 了 晓 
夫 数 Gr 的 增加 ,对流 扰动 的 最 大 振幅 4 的 成 长 率 
变 大 ,线性 成 长 阶段 的 时 间 变 短 。 建 立 了 Pr = 
0.0272 和 Pr = 6. 949 时 , 热 壁面 平均 努 塞 尔 数 Nw 
和 格拉 晓 夫 数 Gr 的 关联 式 。 发 现 随 着 格拉 晓 夫 数 
Gr 数 增加 , 右 壁 面 的 传 热 能 力 增 强 。 


1 数学 物理 模型 


1.1 流体 力学 方程 组 


基于 Boussinesq 假定 ,描述 侧 向 加 热 对 流 问题 
的 流体 力学 方程 组 可 以 表述 为 


VvV:U=0 (1) 

UDU VE+rVU -ay Tg 
ot po 

(2) 

UWT= «VT 0 


式 中 : U(u,0,w) 为 速度 矢量 场 ;7 为 温度 ; 为 压强 ; v 
为 运动 黏 性 系数 ; + 为 时 间 ; « 为 热 扩散 系数 ,x = 了 
和 为 导热 系数 ,C, 为 比 热 ;p 为 密度 ;8 = -ge,,g 为 
重力 加 速度 ;e, 表示 y 轴 单 位 矢量 ,平行 于 y 轴 向 


上 ; a 为 体积 膨胀 系数 ,a = - 全， 下 标 0 
表示 相应 物理 量 的 参考 值 。 
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1.2 边界 条 件 与 初始 条 件 


下 面 给 出 速度 及 温度 的 边界 条 件 。 
当 y=0, 厂 时 ,w=w=0 
当 x=0,d 时 ,wu=w=0 
当 y=0, 厂 时 ,97/9y=0 
当 x=0 时 ,T= 
0 
式 中 :uw 分 别 表 示 水 平流 速 、 垂 向 流速 ;六 为 计算 
区 域 高 度 ;d 为 计算 区 域 宽度 ;7 为 左 壁面 温度 ; 人 
为 右 壁 面 温 度 。 
初始 流速 条 件 为 w = w = 0 , 初始 温度 为 边 壁 
面 突 然 加 温 。 
系统 由 以 下 参数 来 控制 , 即 : 腔 体 的 高 宽 比 4 = 


3 
Tvd ;格拉 晓 夫 数 Gr = SEC AI; 普 朗 特 数 Pr = wx。 
Vv 


为 讨论 方便 ,对 腔 体 尺 寸 进行 量 纲 归 一 化 处 理 , 即 
Y= y/d,X = x/d, 


1.3 ”加 热 条 件 


左 壁 面 温度 保持 恒定 , 右 壁 面 温度 Ts 周期 变 
化 , 即 
Ts 
7 = T+ Ssin( Ty) (4) 


1 
7 = pa 全 十 Thinax ) (5) 


式 中 : AT = Ti -也 
图 1 是 腔 体高 宽 比 4 = 10 时 左右 壁面 的 温度 
分 布 。 


，Tkwa 表示 石壁 面 最 大 温度 ; 


u 
闪 293.15 


图 1 垂直 壁 函 


Fig.1 Temperatures along the vertical walls 


1.4 数值 方法 
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使 用 均匀 交错 网 格 系统 ,网 格 尺寸 分 为 Ax = Ay = 
d/20 和 Ax = Ay = d/30 两 种 。 文 献 [24] 在 Gr = 
5.700 4 x10 情况 下 , 当 Ax = Ay = d/20 时 , 右 壁 面 
努 塞 尔 数 的 平均 值 Nu = 0.0331; 当 Ax = Ay = 
d/30 时 , Nu = 0.033 3 ,得 到 Nu 两 者 的 误差 为 
0.604% ,发 现 特征 量 很 接近 。 因 此 ,采用 Ax = 
Ay =d/20 。 时 间 步 长 为 0.01s。 在 每 个 时 间 步 ,如 
果 节 点 余 量 绝对 值 的 最 大 值 小 于 10”, 并 且 同 一 时 
间 步 内 前 后 两 次 迭代 计算 的 相对 误差 小 于 10“, 认 
为 收敛 。 本 研究 数值 模拟 采用 Fluent 软件 。 


2 侧 向 正弦 加 热 下 对 流 扰动 的 成 长 和 
传 热 特性 


2.1 侧 向 正弦 加 热 下 对 流 扰动 的 成 长 


2.1.1 对 流 扰动 的 成 长 率 

对 于 侧 向 周期 加 热 的 矩形 腔 体 ,突然 加 热 后 观 
察 腔 体 中 对 流 扰动 的 成 长 ,这 时 ,对 流 扰动 的 最 大 振 
幅 4 可 表示 为 

4 cc exp( Ynt) (6) 

式 中 :yy 是 初始 小 扰动 成 长 阶段 的 成 长 率 ; 1 是 扰动 

时 间 。 由 式 (6) 可 以 看 出 , 随 着 扰动 时 间 的 增加 ,对 

流 扰 动 的 最 大 振幅 4,,, 指数 成 长 。 对 流 扰动 的 成 长 
阶段 最 大 振幅 4,, 的 成 长 率 可 被 定义 为 

yn = d(lnA,,)/di (7) 

取 腔 体 中 铅 垂 方向 速度 场 最 大 值 为 对 流 最 大 振 
幅 。 在 对 流 扰动 成 长 的 初始 阶段 最 大 振幅 4,, 的 成 
长 率 应 该 为 常数 。 

2.1.2 普 朗 特 数 Pr = 6. 949 时 对 流 扰动 的 成 长 

对 于 高 宽 比 4 = 10 , 腔 体 宽度 d = 2 cm 的 腔 
体 , 普 朗 特 数 Pr = 6. 949 的 流体 , 侧 向 正弦 周期 加 热 
条 件 下 ,选择 3 个 格拉 晓 夫 数 Gr 从 初始 状态 开始 计 
算 ,观察 腔 体 中 对 流 扰动 的 成 长 。 

图 2 给 出 了 不 同 格拉 晓 夫 数 Gr 情况 下 对 流 扰 
动 成 长 初始 阶段 的 最 大 振幅 4, 随 着 时 间 1 的 变化 。 
在 格拉 晓 夫 数 Gr = 57 004 的 情况 下 ,时 间 1 < 14 s 
时 ,最 大 振幅 4,,, 随 着 时 间 1 线性 变化 。 直 线 的 斜率 
是 一 个 常数 。 也 就 是 说 ,对 流 扰动 的 成 长 阶段 最 大 
振幅 4, 的 成 长 率 是 一 个 常数 。 当 时 间 上 = 14 s 后 ， 


利用 有 限 容积 法 对 流体 力学 方程 组 进行 了 离 
散 ,采用 SIMPLE 算法 求解 流体 方程 组 ,对 流 项 采用 
二 阶 迎 风格 式 , 时 间 项 采用 一 阶 隐 格式 。 计 算 区 域 


最 大 振幅 4,, 随 着 时 间 1 的 变化 偏离 直线 ,最 大 振幅 
4 的 增加 变 小 。 系 统 进入 非 线 性 阶段 。 对 于 不 同 
的 格拉 晓 夫 数 ,直线 的 斜率 是 不 同 的 常数 ,也 就 是 
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说 ,对 流 扰动 的 成 长 阶段 最 大 振幅 4 的 成 长 率 不 
一 样 。 随 着 格拉 晓 夫 数 的 减 小 ,对 流 扰动 的 成 长 阶 
段 最 大 振幅 4,,, 的 成 长 率 变 小 ,线性 成 长 阶段 的 时 
间 变 长 。 通 过 不 同 格拉 晓 夫 数 情 况 下 ,线性 成 长 阶 
段 直线 的 斜率 可 获得 相应 的 成 长 率 。 图 3 给 出 了 成 
长 率 y， 随 格拉 晓 夫 数 Gr 的 变化 。 可 以 看 出 ,对 流 
扰动 的 成 长 阶段 最 大 振幅 4,, 的 成 长 率 随 着 格拉 晓 
夫 数 Gr 的 增加 而 增加 。 对 图 3 的 资料 进行 拟 合 分 
析 ,可 获得 成 长 率 y,, 和 格拉 晓 夫 数 Gr 的 经 验 式 为 
y= 9 x10 Gr , 拟 合 可 决 系数 RR = 0. 963 。 


0.004 
-@- Gr=1.7004 X10’ 
-G0/=2.7004X10" 
0.003 Gr=5.7004X10' 


0.002 


0 6 12 18 24 30 
ts 


图 2 普 朗 特 数 Pr = 6. 949 时 最 大 振幅 随时 间 的 变化 
Fig.2 Variation of the maximum amplitude 


with time at Pr = 6. 949 


1 2 3 4 3 6 
Gr/10’ 


图 3 普 朗 特 数 Pr = 6. 949 时 最 大 振幅 的 成 长 率 
ya 随 格拉 晓 夫 数 Gr 的 变化 
Fig.3 Variation of the growth rate y,, of the 
maximum amplitude with Gr at Pr = 6. 949 
2.1.3 普 朗 特 数 Pr = 0.703 时 对 流 扰动 的 成 长 
对 于 高 宽 比 4 = 10 ,考虑 到 Boussinesd 假定 , 腔 
体 宽度 d = 6 cm 的 矩形 腔 体 , 普 朗 特 数 Pr = 0.703 
的 流体 , 侧 向 正弦 周期 加 热 条 件 下 ,选择 4 个 格拉 晓 
夫 数 Gr 从 初始 状态 突然 加 热 开 始 计算 ,观察 腔 体 中 
对 流 扰动 的 成 长 。 
普天 特 数 Pr = 0. 703 时 不 同 格拉 晓 夫 数 Gr 情 
况 下 对 流 扰动 初始 成 长 阶段 的 最 大 振幅 4 随 着 时 
间 :的 变化 如 图 4 所 示 。 在 格拉 晓 夫 数 Gr = 13 480 
的 情况 下 ,时 间 t < 14 s 时 ,对 流 最 大 振幅 4,, 随 着 
时 间 t 线 性 变化 。 线 性 变化 的 直线 斜率 是 1 个 常数 。 
也 就 是 说 ,对 流 扰动 的 初始 成 长 阶段 成 长 率 是 1 个 
常数 。 当 时 间 上 = 14 s 后 ,最 大 振幅 4,, 随 着 时 间 1 


偏离 直线 变化 ,最 大 振幅 4,, 的 增加 变 小 。 对 流 系 
统 进入 非 线 性 范围 。 当 格拉 晓 夫 数 Gr = 23 480， 
53 480 ,67 398 时 ,对 流 系 统 最 大 振幅 由 线性 转变 为 
非 线 性 成 长 的 时 间 t = 10 s,8 s,6s 。 可 见 , 随 着 格 
拉 晓 夫 数 的 增加 ,对 流 扰 动 的 最 大 振幅 线性 成 长 阶 
段 的 时 间 变 短 。 同 时 可 以 发 现 , 随 着 格拉 晓 夫 数 的 
增加 ,对 流 最 大 振幅 4,, 随 着 时 间 上 直线 变化 的 角度 
变 大 ,直线 的 斜率 变 大 ,也 就 是 说 ,对 流 扰动 的 最 大 
振幅 4,, 的 成 长 率 变 大 。 通 过 不 同 格拉 晓 夫 数 情况 
下 ,线性 成 长 阶段 直线 的 斜率 可 获得 相应 的 成 长 率 。 
线性 成 长 阶段 的 成 长 率 y,, 随 格拉 晓 夫 数 Gr 的 变化 
如 图 5 所 示 。 可 以 看 出 ,对 流 扰动 的 成 长 率 随 着 格 
拉 晓 夫 数 Gr 的 增加 而 增加 。 对 图 5 的 资料 进行 拟 
合 分 析 ,可 获得 对 流 扰动 的 成 长 率 y, 随 着 格拉 晓 夫 
数 Gr 的 变化 关系 式 为 y,, =8 x10 一 cr , 拟 合 可 决 
系数 R? = 0.9077 。 


0.05 
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图 4 普 朗 特 数 Pr = 0.703 时 最 大 振幅 随时 间 的 变化 
Fig.4 Variation of the maximum amplitude 


with time at Pr = 0. 703 


0.3 
全 
有 
A 
0.1 
0 5 4 6 8 
Gr/10’ 
图 5 普 朗 特 数 Pr = 0.703 时 最 大 振幅 的 成 长 率 y， 
随 格拉 晓 夫 数 Gr 的 变化 


Fig.5 Variation of the growth rate y,, of the maximum 


amplitude with Gr at Pr = 0. 703 
2.2 热 壁面 传 热能 力 


正弦 周期 加 热情 况 下 , 普 朗 特 数 Pr = 6. 949 和 
Pr = 0.027 2 时 热 壁 面 平 均 努 塞 尔 数 Nu 的 比较 如 
图 6 所 示 。 可 知 , 对 于 给 定 的 Pr , 随 着 格拉 晓 夫 数 
Gr 的 增 大 , 热 壁面 平均 努 塞 尔 数 Nu 增 大 , 热 壁 面 的 
传 热能 力 越 强 。 对 于 不 同 Pr , Pr 越 大 , 热 壁 面 平均 
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努 塞 尔 数 Nu 增 大 , 热 壁面 的 传 热能 力 越 强 。 对 于 图 
6 中 数据 进行 拟 合 分 析 , 给 出 了 2 种 Pr 下 热 辟 面 平 
均 努 塞 尔 数 Nu 随 着 格拉 晓 夫 数 Gr 变化 的 关系 式 ， 
可 决 系数 及 应 用 条 件 如 表 1 所 示 。 


6 
5 
4 
S 
EE, 
和 Pr-0.027 2 
1 国 忆 /一 0.949 
0 1 和 3 4 5 6 7 
Gr/10’ 
图 6 不 同 Pr 时 的 传 热能 
Fig.6 The heat transfer capacity for different Pr 
表 1 右 壁 面 平均 努 塞 尔 数 随 着 格拉 晓 夫 数 
Gr 变化 的 关系 式 
Tab.1 Correlations for the average Nusselt 
number on the right wall with Gr 
传 热 关联 式 可 决 系数 应 用 条 件 


Nu = 0.00002CrrD R? =0.9936 8.75 x10 < Gr < 6.78 x104 


Nu = 0.00024CrL R? = 0.9284 1.75 x103 < Gr < 1.754 x 10 


3 结 论 


本 研究 通过 对 流体 力学 方程 组 的 数值 模拟 , 研 
究 了 和 矩形 腔 体 侧 向 正弦 周期 加 热 条 件 下 对 流 扰动 的 
成 长 和 传 热 特性 。 可 以 得 到 以 下 结论 。 

1) 对 于 侧 向 正弦 周期 加 热 矩形 腔 体 , 对流 扰动 
的 最 大 振幅 4 的 成 长 率 随 着 格拉 晓 夫 数 cr 的 增 
加 变 大 ,线性 成 长 阶段 的 时 间 随 着 格拉 晓 夫 数 Gr 的 
增加 变 短 。 

2) 对 于 高 宽 比 4 = 10 , 腔 体 宽度 d = 2 cm 的 腔 
体 , 普 朗 特 数 Pr = 6. 949 的 流体 ,对 流 扰动 的 成 长 率 
yw 随 着 格拉 晓 夫 数 Gr 变化 的 经 验 式 为 y。 = 9 x 
10” 6 ;对 于 高 宽 比 4 = 10 、 腔 体 宽度 d = 
6 cm 的 腔 体 ,普度 特 数 Pr = 0.703 的 流体 ,成 长 率 
y 随 着 格拉 晓 夫 数 Gr 的 变化 关系 式 为 yw = 8 x 
10 Gr”, 

3) 当 Pr = 0.0272 时 , 热 壁 面 平均 努 塞 尔 数 Nu 
和 格拉 晓 夫 数 G7 的 关联 式 为 Nu = 0. 000 026r'7 ; 
当 Pr = 6.949 时 , Nu = 0. 000 24Gr'% 。 对 于 不 同 
Pr , Pr 越 大 , 热 壁 面 平均 努 塞 尔 数 Nu 增 大 , 热 壁面 
的 传 热 能 力 越 强 。 
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